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We show that the rate of charged-current events with a hard internal bremsstrahlung photon 
is consistent with the excess candidate events reported by the MiniBoone and LSND (decay in 
flight) experiments. Hard photons radiated by the muon leg in charged-current neutrino interactions 
(i/ M N— >/ryN) are a significant source of background that should be considered by current and future 
Vp-^>v<i neutrino oscillations appearance experiments (e.g. LSND, MiniBoone, SuperK, MINOS, 
T2K and NOVA). 
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Experimental evidence for oscillations among the three 
neutrino generations was reported almost a decade ago. 
The LSND Collaboration [T] has also reported evidence 



for 



oscillations in a A m and mixing angle re- 



gion which is not consistent with the atmospheric and 
solar neutrino oscillation results. This LSND result was 
obtained from the observation of Ve events in a v ^ beam 
(with energies between 20 and 55 MeV) that originates 
from decays of muons at rest [I]. 

The LSND collaboration also searched for v e candi- 
dates in an independent exposure to a beam of i^'s (with 
energies in the range 60 to 200 MeV) from pion decay in 
flight [2]. The observed rate of v e events in the decay 
in flight sample has been reported to be consistent with 
the oscillations parameters from the decay at rest results 
(though with a lower statistical significance). 

The MiniBoone experiment was constructed to inves- 
tigate the origin of the anomalous excess of v e events 
reported by LSND. The and za, beams at MiniBoone 
are at higher energy (0.2-3.0 GeV) 

The first published results 3 from MiniBoone (taken 
with a beam) do not confirm the LSND results. The 
observed yield and energy spectrum of the MiniBoone v e 
candidates (shown in figure 1) is incompatible with the 
oscillation parameters in the LSND allowed region. How- 
ever, the MiniBoone data show an excess of v e candidates 
at energies below their analysis threshold of 475 MeV . 
The total excess which is not understood, is about 150 
events in a sample of « 400,000 charged current 
events, or about 0.03%. When the analysis is extended to 
lower energy jl] (see in figure 2), the excess of v e events 
persists. 

In this communication we show that the calculated 
rate of charged-current neutrino events with a hard 
internal bremsstrahlung photon is consistent with the ex- 
cess v e event candidates reported by MiniBoone. Hard 
photons radiated by the muon leg in charged current 
events ( i/fJSS-^fi'yN) are a significant source of back- 
ground that should be considered by >v e neutrino os- 
cillations appearance experiments. 

The probability for radiation of a hard photon by 
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FIG. 1: The reconstructed neutrino energy (-EJ ec ) distribu- 
tion of the total (top) and excess (bottom) of z/ e charged cur- 
rent quasi-elastic candidate events published by MiniBoone. 
The total excess, which is not understood (below their anal- 
ysis threshold of 475 MeV) is about 150 events in a sample 
of « 400,000 charged current events, or about 0.03% . 
Also shown are the best fit oscillation parameters, and the 
predicted excess from oscillation parameters in the LSND al- 
lowed region. 



muons created in a quasielastic charged current in- 
teraction in the MiniBoone energy range is about 1%. 
The probability to radiate a photon which carries about 
half of the muon energy is about 0.3%. Such high en- 
ergy photons appear as electromagnetic (EM) showers in 
Cerenkov detectors such as MiniBoone. In Cerenkov de- 
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FIG. 3: The the lowest order Born diagram for d 2 o~o/dxdy 
for !^ M N^ j«N 




FIG. 4: The the bremsstrahlung diagram for radiation of a 
hard photon by the muon leg i/^N— > /ryN). 



FIG. 2: The reconstructed neutrino energy (E„ ec ) distribu- 
tion of excess of v e from an updated analysis extended to 
lower energy. 



tectors, muons with energy less than 200 MeV are prac- 
tically invisible, since they are below Cerenkov threhold. 
Therefore, 400 MeV muons, which radiates a 200 MeV 
photon would be invisible in MiniBoone and appear as v e 
candidates. Even if a muon is partially visible in Mini- 
Boone because its energy is above Cerenkov threshold, or 
due to the production of a low light level from the scintil- 
lator in oil (about 25%), its signal would most likely be 
buried within the photon electromagnetic shower. The 
presence of some of these undetected muons may be in- 
ferred in MiniBoone via the observation of a Michel (< 
52.8 MeV) electron from the delayed decay at rest of the 
muon (/z — DAR). However, the efficiency for the detec- 
tion of the presence of such muon events using fj, — DAR 
photons is not 100% for several reasons: 

(1) A fraction of of the DAR electrons are at very low 
electron energy. 

(2) A fraction of the muons decay outside the detector's 
active fiducial volume. 

(3) A fraction of the negatively charged muons are cap- 
tured by atomic nuclei and undergo internal conversion. 

In this communication, we show that the rate of radia- 
tive muon events is significant and should be included in 
the MiniBoone analysis. The details of the inefficiency 
in the detection of invisible or partially visible muons de- 
pend on the specific details of the MiniBoone analysis, 
and are not addressed in this communication. 



The radiative corrections to inclusive charged current 
neutrino scattering has been calculated by by several 
authors [5] including Kiskis, Barlow and Wolfram, and 
Bardin et al. In this communication, we use the lead- 
ing log peaking approximation as derived by Dc Rujula, 
Pctronzio, and Savoy-Navarro [6]. 

The contribution from the muon leg to the double dif- 
ferential cross section in leading log approximation can 
then be written as 
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and d 2 ao /dxdy is the lowest order Bom cross section 
as shown in figure 3. Here, Q 2 — 2ME u xy, M is 
the nucleon mass, y = (E u — E^)/E u , and E^ the en- 
ergy of the outgoing charged muon in the laboratory 
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frame. Eq. ([I]) is expected to be a good approximation 
for \n{2ME v /ml) > ln(l - y) [6]. 

The first term under the integral sign of equation (1) 
comes from the bremsstrahlung diagram of figure 4 
(z/ M N— >/ryN) . Muons of energy £7 M = E^Jz > radi- 
ate photons and contribute to the uncorrected differen- 
tial cross sections for outgoing muons of energy . The 
case of z — > z m corresponds to the case when the initial 
muon originates from quasi-elastic scattering. The case 
of z — ► 1 corresponds to the case where the energy of the 



radiated phone E~ 



0. 



For MiniBoone we make the approximation that the 
cross section is dominated by quasi-elastic scattering 
from nucleons at rest. In that case, the Born cross section 
is a delta function at the quasielastic peak. 



d 2 (T 
dxdy 



G(x,y)S(x-l) 



(3) 



Where G(Q 2 ) is a combination of vector and axial 
form factors. Integrating equation (1) over the elastic 
delta function we obtain the probability dP{y)/dy for a 
muon (from a quasielastic scattering event) with en- 
ergy Efj_- qe to radiate a hard photon and end up in with 
energy E^ = yE v : 



dP{y) 

dy 



a , 2ME v zl l 

^~ ln 2 

Zir mi 



(4) 
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We define r=_E 7 /£' Al _ ge . The case of yi ow =(1 — 
Eu-qe/ 'E v ) is when the energy of the radiated photon 
Ej — ► (r = r low = 0). The case ofy high =\-m[ilE v is 
when the momentum of the final state muon is zero. The 
case of r = r^igh yields the highest possible radiated pho- 



ton energy E 1 — > E 



H-ge 



^ . By Integrating equation 
Q, we obtain the integrated probability P(i? 7 _ m i n ) for 
a muon to radiate a photon with energy _E 7 > £L_ m £„. 

The background for v e events originates from events 
with a high _E 7 . In MiniBoone, the total measured energy 
of an electromagnetic shower for a quasi-elastic events 
with a high photon energy is the sum of the energy 
of the photon and the fraction of energy of the muon 
above Cerenkov threshold. This is given approximately 
by E E m=E 1 + Max (E^ qe - £ 7 -0.2 GeV, 0). 

We now give a few numerical examples. For the case 
of E v = 442 MeV, E^ qe = 414 MeV, and 6> M = 13 de- 
grees (Q 2 =0.01 (GeV/c) 2 ), the probability for radiating 
a hard photon with an energy which is higher than a typ- 
ical experimental muon momentum resolution (e.g. 5%) 
is 0.66%. The probability for radiating a photon with 
E 1 > 207 MeV (i.e. 50% of its energy) is 0.33%. This 
rate an order of magnitude larger than the 0.03% rate of 



excess v e candidates in MiniBoone. Additional numer- 
ical examples at other representative energies are given 
in Table (I) including both lower energy (e.g. 226 MeV 
corresponding to the lower energy range at MiniBoone 
and the upper range of LSND), as well as higher energy 
muon. 

In the MiniBoone experiment the nucleons are bound 
in carbon. For the case of quasi-elastic v scattering, a 
a nucleon of mass Mi ( bound in C 12 ) is scattered to 
a final state nucleon with mass M%. Here, Q 2 — q 2 — 
{E v — Ei) 2 — —mf + 1E v (Ei — pi cos # M ), where E t is the 
final state lepton energy (muon or electron). MiniBoone 
calculates the reconstructed energy (shown in figure 1) 
as follows. 

= 2(M 1 - E B )E l - (E% - 2M 1 E B + m 2 + AM 2 ) 

2[(M 1 -E B )-E l +p l coB9 l ] [ 1 

MiniBoone uses an average removal energy 
Eb = 37MeV. Here, mi is final state lepton mass, pi is 
the final state lepton momentum, and AM 2 = M 2 — M 2 . 
For v e candidates, the final state lepton mass is assumed 
to be the electron mass (which can be neglected) . When 
calculating the reconstructed energy for the radiative 
muon background to the v e candidates, the final state 
energy and momentum of the lepton is Eem Therefore, 
the reconstructed neutrino energy (shown in figures 1 
and 2) is higher than Eem for electromagnetic showers 
at larger angles. 

Although radiative muon events are a large source of 
background in v^^>v e neutrino oscillations appearance 
experiments, radiative corrections have only a small ef- 
fect on the overall quasielastic differential cross section 
at low energies. However, as the precision of the next 
generation neutrino experiments improves, radiative cor- 
rections should be accounted for. It is usual to define 
5 = P (Ey-min = AEfj) as the radiative correction to 
quasielastic scattering. Here A(E I1 ) = a x _E M _ ?e , where 
a is larger than the experimental error on the measure- 
ment of the energy of the final state muon,. The Born 



quasielastic cross section 



dy 



sured quasielastic cross section 
expression. 

d (J quasi d~(7Q—q e 



is related to the mea- 

•3 

L via the following 



dy 



dy 



dy 



1 - S(AE^) 



(6) 



The above expression can be used to correct the value 
of the axial vector mass Ma extracted from quasielastic 
neutrino scattering experiments for radiative effects. In 
general, correcting for radiative effects will increase ex- 
tracted value of Ma (since correcting for radiative effects 
yields is a larger correction at larger Q 2 ). The values of 
the radiative corrections for an experiment with a typ- 
ical muon energy resolution A = 0.05 x E^ qe are also 
given in Table (I). The updated world average value [7] of 
the axial vector mass Ma extracted from quasielastic 
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180.0 


1.051 
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0.442 
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0.714 
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4.573 


13.4 


1.00 


0.0177 


0.0089 


1.703 
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33.3 
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0.0156 


0.0060 


4.001 


4.793 


15.5 


1.40 


0.0182 


0.0088 


1.703 


2.828 


37.6 


2.00 


0.0157 
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0.0085 


10.001 


10.029 


0.6 


0.01 
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0.3 


0.01 


0.0236 


0.0152 


10.001 


10.078 


1.8 


0.10 


0.0201 


0.0125 


20.000 


20.078 


0.9 


0.10 


0.0236 


0.0152 


10.001 


10.133 


2.6 


0.20 


0.0202 


0.0125 


20.000 


20.133 


1.3 


0.20 


0.0237 


0.0152 


10.001 


10.298 


4.0 


0.50 


0.0206 


0.0124 


20.000 


20.297 


2.0 


0.50 


0.0238 


0.0152 


10.001 


10.572 


5.6 


1.00 


0.0210 


0.0123 


20.000 


20.572 


2.8 


1.00 


0.0241 


0.0151 


10.001 


10.791 


6.5 


1.40 


0.0214 


0.0122 


20.000 


20.791 


3.3 


1.40 


0.0243 


0.0150 


10.001 


11.121 


7.7 


2.00 


0.0218 


0.0120 


20.000 


21.120 


4.0 


2.00 


0.0246 


0.0149 



TABLE I: A few numerical examples of the probability for a muon generated quasielastic neutrino scattering to radiate a high 
energy photon with energy greater than 5% and 50% of the energy of the quasielastic muon £ M _ 9e = (pi + ml) 1/2 . The first 
two entries in the table, £ M _ ge = 0.117 GeV (p M = 0.050 GeV/c), and £ M _ 9e = 0.175 GeV (p M =0.140 GeV/c) are in the 
LSND (decay in flight) energy region (and lower range of MiniBoone). The higher energy entries are useful for experiments 
that extract the axial form factor from quasielastic events. 



scattering bubble chamber experiments with deuterium 
is Mf uterium =1.014 ±0.026 GeV/c 2 . For a typical 
energy of 2 GeV, the corrections for radiative effect given 
in Table (I) implies that this average value should be in- 
creased by 8M V j£ « 0.002. For bubble chamber experi- 
ments with antineutrinos, the correction depends on the 
range of Q 2 of the analysis. 

In summary, we show that the background from 
i^N— >/x7N can account for the excess v e candidates in 
the MiniBoone and LSND (decay in flight) experiments. 
Hard photons radiated by the muon leg in charged cur- 
rent events are a significant source of background that 
should be considered by f M — neutrino oscillations ap- 
pearance experiments (e.g. LSND, MiniBoone, SuperK, 
MINOS, T2K and NOVA). 
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